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Abstract. We consider magnetic geodesic flows on the 2-torus. 
We prove that the question of existence of polynomial in momenta 
first integrals on one energy level leads to a Semi-Hamiltonian sys- 
tem of quasi-linear equations, i.e. in the hyperbolic regions the 
system has Riemann invariants and can be written in conservation 
laws form. 



1. Introduction 

In this paper we introduce a new hierarchy of Semi-Hamiltonian sys- 
tem which is naturally related to integrable magnetic flows on surfaces. 

More precisely we consider magnetic geodesic flows on two-torus. 
Consider one energy level and assume it admits a polynomial in mo- 
menta integral of motion. Then we prove that the system of quasi-linear 
equations on the coefficients is in fact Semi-Hamiltonian system. These 
systems were introduced by S.Tsarev and later on studied extensively. 
It is proved in that these systems are integrable by the generalized 
hodograph method. A remarkable theorem by Sevennec [2] states that 
Semi-Hamiltonian property is equivalent to existence of two special co- 
ordinate systems in the space of field variables: Riemann invariants 
and Conservation laws. It is remarkable fact that both forms naturally 
appear for many systems. For example for Benney chains [3j and for 
geodesic flows [1] the Riemann invariants correspond to critical val- 
ues of the integral, and Conservation laws are related to the invariant 
torii of the flow. In this paper the Semi-Hamiltonian property appears 
naturally in the same manner as we shall prove below. 

The problem of existence of integrable systems in the presence of 
gyroscopic forces (which is equivalent to magnetic field) was studied 
by V.V.Kozlov and his students |5],[6],[7]. Topological obstructions to 
the integrability of mechanical systems on surfaces with magnetic fields 
were obtained in [7j. 
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First of all let us recall some facts about geodesic flows on 2-torus 
(without magnetic field). If the geodesic flow is integrable then on the 
torus there are global semi-geodesic coordinates ds"^ = g'^{t, x)dt^ + dx'^ . 
This coordinates can be constructed using invariant Liouville torus 
having the diffeomorphic projection on the base T^. The existence 
of such a torus is proven in [8]. In the semi-geodesic coordinates the 
polynomial in momenta integral has the form 



Pi 



k=0 



-k 



where a 



n-l 



and 



1. The coefficients satisfies the system 



where U 



A 



( 



Vt + A{\J)V^ = 0, 
o^n-i)^) and matrix A has the form 



(1) 



v 


































2a2 — na^ 
3a3 — (n — l)ai 

{n - l)a„_i - 3a„_3 
nan - 2a„_2 



In ^ it is proven that in the case of integrals of degree 3 or 4 in the 
elliptic regions (where matrix A has all eigenvalues different and two of 
them are complex-conjugate) the integral can be reduced to integrals 
of the first or second degree or the metric is fiat. 

We proved in |1] that the system (1) is Semi-Hamiltonian. In this 
paper we generalize this result to the case of nontrivial magnetic field. 



2. The main theorem 

The geodesic fiow on the torus with the Riemannian metric ds^ = 
given by the Hamiltonian equations on T*T 

X = 7^ — , Pj = — t; — r, J = 1,2, 
dpj ox^ 

where H = ^g^-'PiPj- The function F : T*T^ — )■ M is called the first 
integral of the geodesic fiow if 

F = {F,H}, = 0, 

where {F, H}g is Poisson bracket 

^ ' dF dH dH dF 



dx^ dpj dx^ dpj 
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In the case of magnetic geodesic flow the Poisson bracket gets the form 

-vT— — — ^ n(——-—— 

"^^ ~l\(^^^^Pj dx^ dpj J \dp1dp2 dp2dpi 

where fl is the magnetic field. The magnetic geodesic flow given by 
the Hamiltonian equations with respect to magnetic geodesic bracket 

{•; ■}mg 

X' = {x\ H}mg, Pi = {Pi, H}mg- (2) 

Let us choose conformal coordinates (x, y) on the torus. The metric 
becomes the form ds'^ = A{x, y){dx'^ + dy"^). 

We shall consider the problem of existence of polynomial in momenta 

integral of motion F of degree on one energy level H = ^^2a'^ = \- 
We can parameterize the fibres of the energy level over as follows. 
Write 

Pi = \/Acos(^, P2 = sm ip. 
The equations (2) becomes of the form 

1 1 

X — —= cos ip, y — —= sm 

K Ax . ^ 
ip — -p= cosip -p= sm(z? — \ l. 

2A^/A 2A^/A 
The integral F on the energy level has the form 

k=N 
k=-N 

where at = Uk{x, y) + ivk{x, y), a^k = Ofc- 
The condition F = is equivalent to 

{F)x cos (p + {F)y sin (p + F^ (y2A ~ 2X ~ flVA^ = 0. 

We substitute (3) in the last equation and collect terms with respect 
to e''"^. We get 

{ak-i)x + {ak+i)x (Qfc-i)y - («fc+i)j/ A^ i{k - l)afc-i + i{k + l)afc+i 
2 2i 2A 2 

i{k - l)ak-i - t{k + l)ak+i ..^./r- „ ... 
- — tkflV Attk^O, (4) 

where = 0, . . . , + 1 (we assume = at s > TV). 
Let A; be AT + 1 in (4): 

(ajv)a: - N^aN + - ({aN)y - N^aN] = 0. 



2A ^' i V 2A 
Multiply the last identity by A'^/^ 
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(aArA-^/').-i(ajvA-^'/2), = 0. 

Thus oatA"^/^ is a holomorphic function, consequently — A^/^(q;+ 
for some constants a, Taking appropriate rotation in the plane 
{x, y) and dividing F by appropriate constant we can assume that 
a = 1, /3 = 0. Thus we have = A^/^. 

Notice that for = the equation (4) does not contain magnetic 
field and has the form 

(g-i)^ + {ai)x , (a-i)y - {ai)y , K ^'^i ~ ^"-i , Q-i + Qi _ n 
2 + 2l +2A 2 +2A^~""- 
Let us introduce the notation 

^ {0'k-l)x + («A;+l)a: , {0'k-l)y — («fe+l)j/ , 

2 + 2i + 

_^Ay i{k - l)afe_i + i{k + l)ak+i K {k - l)oik-i - (A; + l)ak+i 



2A 2 2A 2 

From equation (4) for k = N we can find the magnetic field: 

n= . 

iA^VAttAf 

Let us substitute this expression of fl into (4) for every k = 1, . . . , N—1. 
We get 

ArgfeA^/2 = A;g^afe. (6) 

The equations (5) and (6) form a system of quasi-linear equations on 
2A'' unknown functions A,uo, ■ ■ ■ , un-i, Vi, . . . , vn-i- This is a quasi- 
linear system of the form 

A{U)U, + B{U)Uy = 0, 

where U = (A, uq, . . . , un^i,vi, . . . , vn-i)'^ ■ We shall write in the last 
section this system explicitly for = 2. Our main result is 

Theorem 1. For any N, the quasi-linear system (5), (6) on coefficients 
of the integral F is Semi-Hamiltonian system. 

3. RiEMANN INVARIANTS 

Consider F as a function on the unite circle -S"^ C C. It is a remark- 
able fact that the critical values of 

F = A^/^z^ + aN-iz""-' + • • • + A^/^^-^ 

are Riemann invariants of the system (5), (6). Indeed, let xi, . . . ,X2n 
be the set of distinct critical points (hyperbolic region) Xi & C C. 
Introduce = F{xk), k — 1, . . . , 2N. From the identity 

F^F^x + Fyv + F^ip^O 
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it follows that Tk are Riemann invariants, because having = we 
are left with F^cos^p + FySin(p — and thus the equation on follows: 



{rk)x + )^k{rk)y = 0, Afe = tantpk, Xk 



e 



Let us check that ri, ...,r2N form a regular coordinate system, that is 
Riemann invariants are functionally independent. Write 

zF' = ^_^^iv-i^. . .+a,z-a.,z-' NA^'/^z 



-N 



Notice that the critical points xi, . . . , X2n are roots of zF' and so by 
Vieta formula 

xi. . .X2N ^ -i- (7) 
For convenience we denote field variables as 

{Hi, . . . , 112n) = (A^/^, aAr_i, . . . , Oo, . . . , ai_Ar). 

Then 

^ ^ + F'ixk)— = —ixk). 

Using (7) we have 

+ l xf-^ ... 



^^M^|^) = (-irdet 



^2N I 1 T,2Ar-l 




Splitting the first column and again using (7) we get 

det [j^^ = (-l)^(xi ...X2N- 1)W = (-1)^+^21^. 

Where W = U^^A the Vandermonde determinant. So we 

have that (/x) <H- (r) is a regular change of variable near every point A 
in the strictly hyperbolic region. 

Remcirk. The field variables A, ao are real, and a_s, s > are 
complex conjugate. Therefore computation of j is equivalent to 

the computation of -TTrrFm — r- 

4. Conservation laws 

The aim of this section is to prove that the system (5), (6) can be 
written in the form of conservation laws. This system has many ex- 
plicit conservation laws. For example, the real part of (4) for k = N 

l—N 

multiplied by A^~ has the form of conservation low 

/ . l-Af \ / . 1-N \ 

Another series of conservation laws can be obtained in the following 
way. The identity (4) at A; = gives us a conservation law 
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Similarly we can get this conservation law for the powers of the integral. 
Namely F"^ generates the conservations law 

where a^"^'* = u^^^ + iv^'^ are corresponding Fourier coefficients of F"^. 
So, we have infinitely many explicit conservation lows. Remarkably 
they are in fact polynomial in the field variables. However we do not 
know if it is possible to get by this method functionally independent 
conservation lows. For this reason we proceed in a different way. We 
show that is possible to generate functionally independent conservation 
laws by invariant tori of the magnetic flow, in a similar way we did it 
in [3] for geodesic flows. 

Imaginary part of (4) for k = N multiplied by A^"^ gives us 

nA = ^ (.„_,A"-""^)^ - ^ (.<«-.A'-"'/^), . (8) 

Let = f{x,y) be a surface invariant under the flow. The invariance 
condition reads as follows: 

fx COS / ^ fy sin / _^ Ax sin / _ A^ cos / ^ _ ^ 

\/A ^/A 2A^/A 2AVA 

or equivalently multiplying by A 

(VAsinf)^- {VAcosf)y + nA = 0. (9) 

Substituting (8) into (9) we get a conservation law corresponding to 
invariant surface ip = f{x,y): 

(^V^sin/ + ^^;^_iA(i-^)/^) - (VA cos / + ^^.^v-iA^^'^)/^) = 0. 

Let us show now that by this method we can get, in the strictly hy- 
perbolic region, 2N conservation laws with functionally independent 

Gk = lm[VAzk + ^ajv-iA(i-^)/2], zk = e'^^K 

Let r be the domain of strict hyperbolicity of F [F' has 2N distinct 
roots on the unite circle). 

Proposition 2. Denote by f the open dense subset ofT defined by the 
condition that ±i is not among the roots of F. Let A* = (/xj, . . . , /i2Ar) 
be any point ofT. Then in a neighborhood of A* there exist 2N func- 
tionally independent conservation laws. 

Corollary 1. Quasi-linear system (5), (6) is Semi-Hamiltonian at any 
point ofT. 
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Proof of Corollary. By Sevennec theorem and the Proposition 1 
the system is Semi-Hamiltonian at any point of F. But it is dense, so 
the condition of Semi-Hamiltonicity extends to the whole F. Applying 
Sevennec theorem again we have that also for A* G F\r there are 
2N independent conservations laws, but wc were not able to construct 
them explicitly. This proves the Corollary. □ 

Proof of Proposition. Define a neighborhood of A* with the help of 
Riemann invariants ri, . . . , r2N- Recall are local coordinates. 

Let rl,k = 1, . . . , 2N are Riemann invariants of A*. That is are 
critical points of F* = {A*)^/^z^+alf_iZ^-^ + - ■ ■+{A*)^/^z-^ . Notice 
that all critical points xl are non-degenerate since we are in the strictly 
Hyperbohc region. We shall assume that odd and even values of the 
index k correspond to maxima and minima respectively. Let e > be 
a small number. Define a neighborhood of F* 

Rs^{F:\n-rl\<^,k^l,...,2N}. 

Here e should be chosen smaller then |min|r^ — r^_,_i|. 

So for any polynomial from it has critical points Xk close to xl 
on the unite circle and critical values close to r^. Define now Zk{s,fj,) 
as solutions of the equations on the unite circle: 

F{zk{e, n), fi) ^ rk + {-!)'' e, 

Zk lies in a neighborhood of Xk- Since we assume that xi,X3, . . . are 
points of maxima and X2,X4, . . . are points of minima then it follows 
that Tfc + (—!)''£ arc regular values for any polynomial F with coef- 
ficients taken from {F should be restricted to a neighborhoods of 
xl ). Important for us is that Zk{s.ji) depend smoothly on e and yu. 
When we shrink e to zero then Zk{fJ^,e) — )■ x^ by the construction. In 
the following we shall decrease e, but keeping away from zero, in order 
to get the needed neighborhood R^. 
Define the following functions on R^: 

Gk = lm[^Zk + -^a^_iA(i-^)/2] ^ 



2N 



— Im 
2i 



(where Zk depend on iis implicitly). We claim that one can choose 
£ > sufficiently small in order to have 

d(//i, . . . ,/U2iv) 

This would imply the claim. We have 

dF I 



dCk I _ VA* / 1 \ 9^,1 iM=^' , 
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Where Rli contains all terms of explicit derivation of 
We have 



det 



/dG 



2N 



2i 



xdet 



dm 



Mention that when e — >■ 0, 




1+ 



which is finite since Xk i 



by assumptions. Also zl{e) — )■ xl and so F'{ii*, z^) 0. Therefore the 
determinant in brackets when £ — )■ tends to the determinant of the 
matrix 

'dFXxU^Y 

dm 

This is exactly the determinant of the matrix considered in the section 
3. So it is equal —2W{xl^ . . . ^x^j^) and does not vanish. But then it 

follows that for small e > det ^|^(^)^ 7^ in a neighborhood of A*. 
Proposition 2 and Theorem 1 are proved. 

5. Discussion 
In this section we discus some open problems. 

In the case of n = 2 the equations (5), (6) on functions A,uo,Ui,Vi 
have the form 

[VAui)^ - i^/Avi)y = 0, 




2VA 

Introduce the new functions /, g 

f 



0, 



0. 



fx + 9y = 0, 

(/A), - {gA)y = 0, 
{uo)x + 2A^ - -g{fy - g^) = 0, 



We have 



-{uo)y + 2A^ + -f{fy - g^) = 0. 
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0, 



A 



B 



This system can be written in the form 

A{U)U, + B{U)Uy 

where U = {A,uoJ,g), 

/ 1 \ 

/OA 
2 1 

\o -i/y 

The problem of existence of periodic solution is very interesting. The 
systems of such form (non-evolution form) were considered in (TU] from 
the point of view of blow-up analysis along characteristic curves. It 
would be very interesting to apply these ideas to our system. 
Notice 

gfy = {f9)y - f9y = {f9)y + ffx = 0, 
f9x = if9)x - 9 fx = {f9)x + 99y 



1 











1 


\ 




-9 








-A 





















2 


-1 


-f 





/ 



So we have 



(uo)x + 2A^ + ^{gg^ - //^) - ^{fg)y 



0. 



Thus we have explicit conservation laws form for our system 

fx + 9y = 0, 
(/A),. - (^A), = 0, 



{u, + 2K+-{g^-f))^ 



,1 



f9)y = 0, 



[-u, + 2A-\{g'-f))y-{\fg)^ 



0. 



(10) 



Lust two conservation lows are very interesting. Let us recall that a 
hyperbolic diagonal system 



{ri)x + Ai(ri, 
is Semi-Hamiltonian if 



'-^j'- ki 



fn){ri)y = 0,i = l, ..,n 



where 



-pk 

ki 



It can be proved that a diagonal system is Semi-Hamiltonian if and only 
if it can be written in some coordinates as a system of conservation laws 
(see [2]). We have the diagonal metric (see P) gu = Hf, and Lame 
coefficients can be found from the over-determined system 



dk In Hi 



ik • 
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By Pavlov-Tsarev theorem [TT] if the Semi-Hamiltonian system has 
two conservation lows of the form 

Fx + Gy = 0, Fy + = 0, 

then the corresponding metric gu is Egorov metric i.e. 

In such the metric is potential: gu = dia{r) for a function a. So it 

follows that our system for = 2 is Egorov Semi-Hamiltonian system, 
since we have two conservation lows (10), (11). By similar calculations 
one can check that for = 3 our system is also Egorov system. It 
would be interesting to prove this fact for arbitrary A^. 

Another interesting problem is to find Poisson bracket of hydrody- 
namical type for the system in the form of Dubrovin-Novikov [12j or 
Ferapontov-Mokhov [T3] . 
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